Stacking-order dependence in thermoelectric transport of biased trilayer graphene 
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We numerically study the thermoelectric and thermal transport in trilayer graphene with different 
stacking orders in the presence of interlayer bias under a strong perpendicular magnetic field. In 
biased ABA-stacked case, we find that the thermoelectric conductivity displays different asymptotic 
behaviors with the varying of the temperature, similar to that of monolayer graphene. In the 
high temperature regime, the transverse thermoelectric conductivity a xy saturates to a universal 
value 2.77kBe/h at the center of each LL, while it displays a linear temperature dependence at low 
temperatures limit. The calculated transverse thermal conductivity K X y exhibits two plateaus away 
from the band center. The transition between the two plateaus is continuous, which is accompanied 
by a pronounced peak in the longitudinal thermal conductivity k xx - In biased ABC-stacked case, it is 
found that both the thermoelectric conductivity and thermal conductivity have similar properties to 
the biased bilayer graphene, which is consistent with the behavior of a band insulator. The obtained 
results demonstrate the sensitivity of the thermoelectric conductivity to the band gap near the Dirac 
point. We also verify the validity of the Mott-relation and the generalized Wiedemann-Franz law. 

PACS numbers: 72.80.Vp; 72.10.-d; 73.50.Lw, 73.43.Cd 



I. INTRODUCTION 

Recently, much attention has been paid to the ther- 
moelectric transport properties of graphene both ex- 
perimentallj*i~— and theoretically 4 " 8 -. In experiments, 
the thermopower (the longitudinal thermoelectric re- 
sponse) and the Nernst signal (the transverse response) 
in the presence of a strong magnetic field are found 
to be large, reaching the order of the quantum limit 
fcg/e, where ks and e are the Boltzmann constant and 
the electron charge, respectivelyi~— . Besides monolayer 
graphene, bilayer graphene is also very interesting. Ex- 
periments have shown that bilayer graphene exhibits 
tunable bandgap2*i£ in the presence of an applied bias 
voltage, exhibiting similar properties to semiconductors. 
Thermoelectric measurement^ shows that the room- 
temperature thermopower with a bias voltage can be en- 
hanced by a factor of 4 compared to that of the monolayer 
graphene or unbiased bilayer graphene, making it a more 
promising candidate for future thermoelectric applica- 
tions. Theoretical calculations from the tight-binding 
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models for monolayer and bilayer grapheni 
agreement with the experimental observations^ - — i 11 t 14 . 

More recently, the focus of the study of graphene sys- 
tems has gradually extended to trilayer graphene^ - —. In 
trilayer graphene, the interlayer coupling in stacked lay- 
ers of graphene gives rise to even richer electronic trans- 
port properties. Experimental and theoretical studies 
show thati^ 3 - - — , the electronic structure and the Landau 
level (LL) spectrum at the vicinity of the Dirac point 
are very sensitive to the stacking order of the graphene 
layers. Trilayer graphene has two stable stacking orders: 
(1) ABA (Bernal) stacking, where the atoms of the top- 
most layer lie exactly on top of those of the bottom layer; 



and (2) ABC (rhombohedral) stacking, where atoms of 
one sublattice of the top layer lie above the center of 
the hexagons in the bottom layer. This seemingly small 
distinction in stacking order results in a dramatic dif- 
ference in band structures. The low-energy band struc- 
ture for ABA-stacked trilayer graphene contains both lin- 
ear and hyperbolic bands, similar to the combined spec- 
trum of monolayer graphene and bilayer graphen e 27 ' 28 , 
while ABC-stacked case presents approximately cubic 
bands 25 . Moreover, the LL spectrum of ABA-stacked 
case in a perpendicular magnetic field B can be viewed 
as a superposition of VS-dcpendent monolayer-like LLs 
and B-dependent bilayer-like LL o 23 ' 29 ' 30 . On the other 
hand, the LLs of ABC -stacked case are given by E n oc 
B z l 2 ^Jn(n — l)(n — 2) with Berry's phase 3ir£t-&. In- 
terestingly, when a bias voltage or a potential difference 
is applied to the top and bottom graphene layers, ABA- 
stacked case exhibits a semi-metallic band structure with 
a tunable band overlap between the conduction and va- 
lence band a 31 ' 32 , whereas ABC-stacked case exhibits a 
semiconducting band structure with a tunable band gap, 
similarly to bilayer graphen o 23 ' 25 ' 26 ' 33 . Owing to their 
distinctive band structures, ABA- and ABC-stacked tri- 
layer graphene are expected to exhibit rich novel thermo- 
electric transport properties. However, theoretical un- 
derstanding of the thermoelectric transport properties of 
trilayer graphene is limited compared to that of mono- 
layer or bilayer graphene. In particular, the influence of 
different stacking orders on the thermoelectric transport 
properties has not been studied so far, which is highly 
desired. 

In this paper, we carry out a numerical study of the 
thermoelectric transport properties in both ABA- and 
ABC-stacked trilayer graphene systems in the presence of 
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electrostatic bias between the top and bottom graphene 
layers. The effects of disorder and thermal activation 
on the broadening of LLs are considered. In biased 
ABA-stacked case, the thermoelectric coefficients exhibit 
unique characteristics near the central LL due to the LL 
crossing of electron and hole bands, which are quite dif- 
ferent from those of biased bilayer graphene. Both the 
longitudinal and the transverse thermoelectric conduc- 
tivities are universal functions of the effective bandwidth 
and temperature, and display different asymptotic be- 
haviors in different temperature regimes. The Nernst 
signal displays a peak at the central LL with a height 
of the order of fcs/e, and changes sign near other LLs, 
while the thermopower behaves in an opposite manner. 
The peak values of the Nernst signal and thermopower 
are very large, compared with monolayer graphene due to 
the semi-metallic band overlap near zero energy. The va- 
lidity of the semiclassical Mott relation is found to remain 
valid at low temperatures. In biased ABC-stacked case, 
we observe quite different behavior from biased ABA- 
stacked case near the central LL. Around the Dirac point, 
the transverse thermoelectric conductivity exhibits a pro- 
nounced valley at low temperatures. This is attributed to 
the opening of a sizable gap between the valence and con- 
duction bands in biased ABC-stacked case. In addition, 
we have calculated the thermal transport coefficients of 
electrons for both biased ABA- and ABC-stacked trilayer 
graphene systems. In biased ABA-stacked case, the cal- 
culated transverse thermal conductivity n xy exhibits two 
plateaus away from the band center. The transition be- 
tween the two plateaus is continuous, which is accom- 
panied by a pronounced peak in the longitudinal ther- 
mal conductivity n xx . In biased ABC-stacked case, the 
transverse thermal conductivity n xy displays an appar- 
ent plateau with n xy — 0, which is accompanied by a 
valley which provides an additional evidence for 

the band insulator behavior. We further compare the cal- 
culated thermal conductivities with those deduced from 
the Wiedemann- Franz law, to check the validity of this 
fundamental relation in trilayer graphene systems. 

This paper is organized as follows. In Sec. II, we in- 
troduce the model Hamiltonian. In Sec. Ill and Sec. 
IV, numerical results based on exact diagonalization and 
thermoelectric transport calculations are presented for 
biased ABA-stacked and ABC-stacked trilayer graphene 
systems, respectively. In Sec. V, numerical results for 
thermal transport coefficients are presented. The final 
section contains a summary. 



II. MODEL AND METHODS 

We consider a trilayer graphene system consisting of 
three coupled hexagonal lattices including inequivalent 
sublattices A\, B\ on the bottom layer, A2, -B2 on 
the middle layer, and A 3 , B 3 on the top layer. The 
three graphene layers are arranged in the ABA (Bernal) 
or ABC (rhombohedral) stacking orders^, as shown in 



(a)ABA (b)ABC 




FIG. 1: (color online). Schematic of trilayer graphene lattice 
with ABA and ABC stacking, where the blue/red/green lines 
indicate links on the bottom/middle/top layers. 

FigJT] The difference between ABA and ABC stacking 
is the top layer. For ABA stacking, the top layer will 
be exactly above the bottom layer without any relative 
shift. For ABC stacking, A3 sublattice of the top layer 
lies above the centers of the hexagons in the bottom layer, 
S3 sublattice lies above the centers of the hexagons in 
the middle layer. Here, the in-plane nearest-neighbor 
hopping integral between Ai and Bi atoms is denoted 
by jAiBi = 7o with i = l,---,3. For the interlayer 
coupling, we take into account two largest hopping in- 
tegrals. For ABA stacking, the largest interlayer hop- 
ping is between a B\ (B3) atom and the nearest A2 
atom r )BiA 2 = IB3A2 — 7i ■ The smaller hopping is 
between a B2 atom and three nearest A± (A3) atoms 
1b 2 Ax — 1b 2 a 3 — 73- For ABC stacking, the largest 
interlayer hopping is between a B± (B2) atom and the 
nearest A2 (A3) atom JbiA? — 1b 2 a 3 — 71 ■ The smaller 
hopping is between a B2 (B3) atom and three nearest 
Ai (A2) atoms ^b 2 Ai — 1b 3 a 2 — 73- The values of these 
hopping integrals are taken to be 70 = 3.16 eV, 71 = 0.39 
eV, and 73 = 0.315 eV^. 

We assume that each monolayer graphene has totally 
L y zigzag chains with L x atomic sites on each zigzag 
chain 36 .. The size of the sample will be denoted as 
N = L L z , where L z — 3 is the number of 

monolayer graphene planes along the z direction. Wc 
have confirmed that the calculated results does not de- 
pend on the system sizes (as long as the system lengths 
are reasonably large}^. In the presence of an applied 
magnetic field perpendicular to the plane of biased tri- 
layer graphene, the lattice model for ABA stacking in 
real space can be written in the tight-binding form: 

H = -70 (Y, eiaWc L c UV + eiaiic 2i„C2j° 

(ij)cr {ij)cr 
(ij)io- 

+ J2 ela * JC h„B 3 C2i<?A 2 )-J 3 {J2 ^^IjaB^UaA, 

+ emi 4j<TB 2 C 3i<rA 3 ) + h.C. + ^ Wi(c\ ia C U(T 

(ij) 3 a i<? 
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i(T c 2ia + 4 iCT c 3 i CT ) + ^2(eic\ ia c lia + e 2 4i„C3iatt) aji (E F ,T) = — f 
%a eT J 



where c 



(4 



4»<mJ> 4nja ( c lnja B J are creation op- 
erators on A m and B m sublattices in the m — th layer 
(to = 1, • ■ ■ , 3), with a as a spin index. The sum YlujSa 
denotes the intralayer nearest-neighbor hopping in three 
layers, Yluj) lC r stands for the interlayer hopping between 
the Bi (B 3 ) sublattice in the bottom (top) layer and the 



A 2 sublattice in the middle layer, and Yltij} 3 a stands 
for the interlayer hopping between the B 2 sublattice in 
the middle layer and the A\ (A3) sublattice in the bot- 
tom (top) layer, as described above. For the biased sys- 
tem, the top and the bottom graphene layers gain differ- 
ent electrostatic potentials, and the corresponding energy 
difference is given by A g = e 2 — ei where e\ = —^A g , 
and e 2 = \A g . For illustrative purpose, a relatively large 
asymmetric gap A g = O.370 is assumed, w; is a random 
disorder potential uniformly distributed in the interval 
Wi G [—W/2,W/2] r yo. The magnetic flux per hexagon 



= E 



o 



= j^- is proportional to the strength of the 
applied magnetic field B, where M is assumed to be an 
integer and the lattice constant is taken to be unity. 

For ABC-stacked trilayer graphene in the presence of 
bias voltage, the Hamiltonian can be written as: 



H = -7o(£ eW 'M- c ^ 



(ij)a- 

(ij)ia (ij)so- 

+ Y e ' ta " J C 3jaB 3 C ^A 2 ) + h.C. + Y Wi{c\ l(T C lia 
(13)30 ia 

+ c\ ia c 2l(T + 4 ia c 3lry ) + ^{tic\ ia cuv + e 2 cl ia c 3la l2) 

ia 

The sum J2(ij)a denotes the intralayer nearest-neighbor 
hopping in three layers, Y^Uj)^ stands for the interlayer 
hopping between the B\ (B 2 ) sublattice in the bottom 
(middle) layer and the A 2 (A3) sublattice in the middle 
(top) layer, and Yl(ij) 3 a stands for the interlayer hopping 
between the B 2 (B 3 ) sublattice in the middle (top) layer 
and the A\ (A 2 ) sublattice in the bottom (middle) layer, 
as described above. 

In the linear response regime, the charge current in 
response to an electric field and a temperature gradient 
can be written as J — <tE + &(— VT), where a and a 
are the electrical and thermoelectric conductivity tensors, 
respectively. The transport coefficient be calcu- 

lated by Kubo formula and a xx can be obtained based 
on the calculation of the Thouless number—. In practice, 
we first calculate the T = conductivities aji(E F ), and 
then use the relational 



aji(E F ,T) = J deaji(e 



dm 

de 



to obtain the finite-temperature electrical and ther- 
moelectric conductivity tensors. Here, f(x) = 
l^ e (x-E F )/k B T _|_ j s tjjg p erm i distribution function. 

At low temperatures, the second equation can be ap- 
proximated as 



aji(E F ,T) 



ir 2 k 2 B T da Jt (e,T) 



3e 



dc 



(4) 



E = Ef 



which is the semiclassical Mott relatio n 37 i 38 . The ther- 
mopower and Nernst signal can be calculated subse- 
quently from 3 ^ 



S X y 



E x 
E v 



— Pxx^xx Pyx™yxi 



PxxQyx H~ PyxOt-xx- 



(5) 



The thermal conductivity, measuring the magnitude 
of the thermal currents in response to an applied tem- 
perature gradient, which usually includes electron and 
phonon contributions. In our numerical calculations, 
phonon-rclatcd thermal conductivity is omitted. The 
electronic thermal conductivities Kji at finite tempera- 
ture assume the forms 3 ^ 



dm 

de 



K Jt (E F ,T) = -^fj dea ]t (e)(e^E F ) 2 (^ 

- Tciji(E F , T)<jJ i 1 (E F , T)ciji(E F , T). (6) 

For diffusive electronic transport in metals, it is estab- 
lished that the Wiedemann- Franz law is satisfied between 
the electrical conductivity a and the thermal conductiv- 
ity k of electrons^: 



oT 



L. 



(7) 



where L is the Lorentz number and takes a constant 
value: L = ^(^f) 2 - The validity of this relation will 
be examined for the present trilayer graphene. 



III. THERMOELECTRIC TRANSPORT IN 
BIASED ABA-STACKED TRILAYER 
GRAPHENE SYSTEMS 

We first show the calculated thermoelectric conductiv- 
ities at finite temperatures for biased ABA-stacked tri- 
layer graphene. As shown in Figl^a) and (b), the trans- 
verse thermoelectric conductivity a xy displays a series of 
peaks, while the longitudinal thermoelectric conductiv- 
ity a X x oscillates and changes sign at the center of each 
LL. At low temperatures, the peak of a xy at the central 
LL is higher and narrower than others, which indicates 
that the impurity scattering has less effect on the central 
LL. These results are qualitatively similar to those found 
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FIG. 2: (color online). Thermoelectric conductivities at finite temperatures of biased ABA-stacked trilayer graphene. (a)- 
(b) a xy (Ep,T) and oi xx (Ef,T) as functions of the Fermi energy Ef at different temperature T. (c) shows the temperature 
dependence of a xy (EF,T) for trilayer graphene. (d) Comparison of the results from numerical calculations and from the 
generalized Mott relation at two characteristic temperatures, UbT/Wl = 0.1 and UbT/Wl = 1.5. Here, the width of the 
central LL Wl/^o = 0.0069. The asymmetric gap A 9 = O.370. The system size is taken to be N — 96 x 24 x 3, the magnetic 
flux <f) — 2n/48, and the disorder strength w = 0.1. 



in monolayer graphene, but some differences exist due to 
the overlap of electron and hole bands. At low temper- 
atures, more oscillations are observed in the higher LLs 
than monolayer graphene, in consistent with the further 
lifting of the LL degeneracy in biased ABA-stacked case. 
As shown in Fig[2Jb) , around zero energy, the peak value 
of a xx shows different trends with increasing tempera- 
ture. It first increases with T in the low-temperature 
region, and then decreases with T at high tempera- 
tures. This is due to the competition between - T and 

rf<Tj 'ie ' T ' > °f Eq.Q. The peak value of a xx could either in- 
crease or decrease depending on the relative magnitudes 
of these two terms. At high temperatures, <7jj(e,T) be- 
comes smooth, and consequently a xx begins to decrease. 
In FigfSJc), we find that a xy shows different behavior 
depending on the relative strength of the temperature 
ksT and the width of the central LL Wl (Wl is deter- 
mined by the full-width at the half-maximum of the a xx 
peak). When k B T -C W L and Ep -C W L , a xy shows lin- 
ear temperature dependence, indicating that there is a 
small energy range where extended states dominate, and 
the transport falls into the semi-classical Drude-Zener 
regime. When Ep is shifted away from the Dirac point, 



the low-energy electron excitation is gapped due to An- 
derson localization. When k B T becomes comparable to 
or greater than Wl, a xy for all LLs saturates to a con- 
stant value 2.77k B e/h. This matches exactly the uni- 
versal value (\ri2)k B e/h predicted for the conventional 
integer quantum Hall effect (IQHE) systems in the case 
where thermal activation dominate s 37 ' 38 , with an addi- 
tional degeneracy factor 4. The saturated value of a xy in 
biased ABA-stacked case is in accordance with the four- 
fold degeneracy at zero energy. In the presence of bias 
voltage, the valley degeneracy of the LLs usually is lifted 
by the interlayer potential asymmetry, so that the 12-fold 
energy levels (four and eight levels from the monolayer- 
like and the bilayer-like subbands, respectively) split into 
six different levels with twofold spin degeneracy. How- 
ever, near the Dirac point, the interlayer potential asym- 
metry causes hybridization of the linear and parabolic 
chiral bands, which leads to the fourfold degeneracy for 
zero energy Landau levels ^2r— 

To examine the validity of the semiclassical Mott rela- 
tion, we compare the above results with those calculated 
from Eq.(TJ}, as shown in Figl^fd). The Mott relation is 
a low-temperature approximation and predicts that the 
thermoelectric conductivities have linear temperature de- 
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FIG. 3: (color online). Thermoelectric conductivities at finite temperatures of biased ABC-stacked trilayer graphene. (a)-(b) 
a xy (EF, T) and a X x(Ep, T) as functions of Fermi energy at different temperatures, (c) Comparison of the results from numerical 
calculations and from the generalized Mott relation at two characteristic temperatures, fcsT/A 9 = 0.001 and ksT/Ag = 0.05. 
The asymmetric gap A 9 = O.370. The other parameters are chosen to be the same as in Fig. 2. 



pendence. This is in agreement with our low-temperature 
results, which proves that the semiclassical Mott relation 
is asymptotically valid in the Landau-quantized systems, 
as suggested in Ref. 



IV. THERMOELECTRIC TRANSPORT IN 
BIASED ABC-STACKED TRILAYER 
GRAPHENE SYSTEMS 

For biased ABC-stacked trilayer graphene, we show the 
calculated a xx and a xy at finite temperatures in FigJ3] 
As seen from Fig(3Ja), ct xy displays a pronounced valley 
at low temperature, in striking contrast to ABA-stacked 
case with a peak at Ep = 0. These results are qualita- 
tively similar to those found in biased bilayer graphene^. 
This behavior can be understood as due to the split of 
the valley degeneracy in the central LL by an opposite 
voltage bias added to the top layer and the bottom layer. 
This is in consistent with the opening of a sizable gap be- 
tween the valence and conduction bands in biased ABC- 
stacked trilayer graphene^. a xx oscillates and changes 
sign around the center of each split LL. In FigJ^c), we 
also compare the above results with those calculated from 
the semiclassical Mott relation using Eq. (£[]). The Mott 



relation is found to remain valid at low temperatures. 

We further calculate the thermopower S xx and the 
Nernst signal S xy . In Fig|4ja)-(b), we show the cal- 
culated S xx and S xy in biased ABA-stacked trilayer 
graphene. As we can see, S xy (S xx ) has a peak (peaks) 
at the central LL (the other LLs), and changes sign near 
the other LLs (the central LL). At zero energy, both p xy 
and a xx vanish, leading to a vanishing S xx . Around zero 
energy, because p xx a xx and p xy ot xy have opposite signs, 
depending on their relative magnitudes, S xx could either 
increases or decreases when Ep is increased passing the 
Dirac point. In biased ABA-stacked CclSG, Sxx is domi- 
nated by p xy a xy , consequently, S xx increases to positive 
value as Ep passing zero. At low temperatures, the peak 
value of S xx near zero energy is ±0.81£;_B/e (±69.8 /W/K) 
at ksT — OAWl, which is in agreement with the mea- 
sured value (±70 /1V/K). 42 With the increase of temper- 
ature, the peak height increases to ±3.75ks /e (±323.14 
/iV/K) at k B T = 1.5W L . On the other hand, S xy 
has a strong peak structure around zero energy, which 
is dominated by p xx a xy . The peak height is 7.82fcs/e 
(673.85 fiV/K) at k B T = 1.5W L . The large magnitude 
of S xy and S xx near zero energy can be attributed to 
the semi-metal type dispersion of biased ABA-stacked 
trilayer graphene, and the fact that the system is in the 
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FIG. 4: (color online). The thermopower S X x and the Nernst signal S xy as functions of the Fermi energy in (a)-(b) biased 
ABA-stacked trilayer graphene, and (c)-(d) biased ABC-stacked trilayer graphene at different temperatures. The parameters 
in these two systems are chosen to be the same as in Fig. 2 and Fig. 3, respectively. 



vicinity of a quantum Hall liquid to insulator transition, 
where the imbalance between the particle and hole types 
of carriers should be significant. The thermoelectric ef- 
fects are very sensitive to such an imbalance in Dirac 
materials in comparison with conventional metals. 



In FigJUc)-(d), we show the calculated S xx and S xy in 
biased ABC-stacked case. As we can see, S xy (S xx ) has 
a peak (peaks) around zero energy (the other LLs), and 
changes sign near the other LLs (zero energy). These 
results are qualitatively similar to those found in biased 
ABA-stacked case. In our calculation, we find that S xx 
is always dominated by p xx a xx , consequently, S xx de- 
creases to negative values as Ep passing zero. This is 
different from biased ABA-stacked case. At low tem- 
peratures, the peak value of S xx near zero energy is 
±0.86fc B /e (±74.11 fiV/K) at k B T = 0.01 A g . With 
the increase of temperature, the peak height increases 
to ±2.23fc B /e (±192.16 fiV/K) at k B T = 0.1A 9 . On the 
other hand, S xy has a peak structure around zero en- 
ergy, which is dominated by p xy oi xx . The peak height is 
4.69fc s /e (404.14 fiV/K) at k B T = 0.lA g . 



V. THERMAL CONDUCTIVITY IN BIASED 
ABA- AND ABC-STACKED TRILAYER 
GRAPHENE SYSTEMS 

We now focus on thermal conductivities. In Fig. [SJ 
we show results of the transverse thermal conductivity 
n xy and the longitudinal thermal conductivity k xx for bi- 
ased ABA-stacked trilayer graphene at different temper- 
atures. As seen from FigJSfa) and (b), n xy exhibits two 
flat plateaus away from the central LL. The values of the 
plateaus ±0.00457 fc s /h (±0.048 nW/(K-m)) 

at k B T = O.IWl- With the increase of temperature, the 
values of the plateaus increase to ±0.0237o£;b//i (±0.24 
nW/(K-m)) at k B T = 0.5VLl- At low temperatures, 
the transition between these two plateaus is smooth and 
monotonic, while at higher temperatures, K xy exhibits 
an oscillatory feature at k B T — 0.5Wl between two 
plateaus. On the other hand, k xx displays a peak near 
the center LL, and its peak value increases quickly with 
T. The peak height is 0.0287 fci3//i (0.3 nW/(K-m)) at 
k B T = 0.5VLl. To test the validity of the Wiedemann- 
Franz Law, we compare the above results with those 
calculated from Eq.([7]), as shown in FigfSJc) and (d). 
The Wiedemann-Franz Law predicts that the ratio of 
the thermal conductivity k to the electrical conductiv- 
ity a of a metal is proportional to the temperature. This 
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FIG. 5: (color online), (a)-(b) Thermal conductivities k, xv {Ef ,T) and k X x(Ef,T) as functions of the Fermi energy in biased 
ABA-stacked trilayer graphene at different temperatures, (c)-(d) Comparison between the thermal conductivity as functions 
of the Fermi energy from numerical calculations and from the Wiedemann- Franz Law at two characteristic temperatures. The 
parameters used here are the same as in Fig. 2. 



is in agreement with our low-temperature results, but 
apparent deviation is seen at higher temperatures. 

In Fig. [SI we show the calculated thermal conductivi- 
ties k xx and K xy for biased ABC-stacked case. As seen 
from Figj6ja) and (b), around zero energy, a flat re- 
gion with n xy = is found at low temperatures, which 
is accompanied by a valley in k xx . These features are 
clearly in contrast to those of ABA-stacked case due 
to the presence of an energy gap between the valence 
and conduction bands. When temperature increases to 
ksT = 0.06A g , the plateau with n xy = disappears, 
while k xx displays a large peak. The peak height k xx 
is 0.Uj k B /h (1.49 nW/(K-m)) at k B T = 0.06A g . In 
FigMc) and (d), we also compare the above results with 
those calculated from the Wiedemann-Franz Law using 
Eq.([7]). We find that the Wiedemann-Franz Law remain 
valid at low temperatures. 

VI. SUMMARY 

In summary, we have numerically investigated the ther- 
moelectric and thermal transport properties of biased tri- 
layer graphene with different stacking orders in the pres- 
ence of both disorder and a strong magnetic field. In 
biased ABA-stacked case, the thermoelectric coefficients 



exhibit unique characteristics due to the LL crossing of 
electron and hole bands that are strongly suggestive of a 
semi-metallic band overlap. We find that the thermoelec- 
tric conductivities display different asymptotic behavior 
depending on the ratio between the temperature and the 
width of the disorder-broadened LLs, similar to those 
found in monolayer graphene. In the high temperature 
regime, the transverse thermoelectric conductivity a xy 
saturates to a universal value 2.77fcse/ft. at the center of 
each LL, and displays a linear temperature dependence 
at low temperatures. The calculated Nernst signal S xy 
shows a strong peak at the central LL with heights of 
the order of fcs/e, and changes sign at the other LLs, 
while the thermopower S xx has an opposite behavior. 
The calculated transverse thermal conductivity n xy ex- 
hibits two plateaus away from the band center. The tran- 
sition between these two plateaus is continuous, which is 
accompanied by a pronounced peak in longitudinal ther- 
mal conductivity k xx . The validity of the Wiedemann- 
Franz law relating the thermal conductivity k and the 
electrical conductivity a is verified to be valid only at 
very low temperatures. 

In biased ABC-stacked case, the thermoelectric co- 
efficients display quite distinct behaviors from those of 
ABA-stacked case. Around the Dirac point, the trans- 
verse thermoelectric conductivity a xy exhibits a pro- 
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FIG. 6: (color online), (a)-(b) Thermal conductivities k xv (Ef,T) and k xx (Ef,T) as functions of the Fermi energy in biased 
ABC-stacked trilayer graphene at different temperatures, (c)-(d) Comparison between the thermal conductivity as functions of 
the Fermi energy obtained from numerical calculations and from the Wiedemann- Franz Law at two characteristic temperatures. 
The parameters chosen here are the same as in Fig.3. 



nounced valley with a xy = at low temperatures, in 
striking contrast to ABA-stacked case with a peak. The 
validity of the semiclassical Mott relation between the 
thermoelectric and electrical transport coefficients is ver- 
ified to be satisfied only at very low temperatures. Fur- 
thermore, the transverse thermal conductivity n xy has a 
pronounced plateau with n xy = 0, which is accompanied 
by a valley These are consistent with the open- 

ing of sizable gap between the valence and conductance 
bands in biased ABC-stacked case. 

We mention that in our numerical calculations, the 
flux 2ir/M in each hexagon gives a magnetic field of the 
strength B ~ 1.3 X 10 5 /M Tesla^. Thus the magnetic 
field B we used is about 2700 Tesla. This magnetic field 
is much stronger than the ones which can be achieved in 
the experimental situation, as limited by current compu- 
tational capability. In our calculation, the system size is 
taken to be N = 96x24x3, and M is taken to be L x or L y 



in consistence with periodic boundary conditions, which 
limits us to extremely strong magnetic fields. However, 
the obtained thermoelectric transport coefficients exhibit 
universal behaviors, as long as M is not too small (greater 
than 10). 
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